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Magnetic configurations in cubic Bi2MnFe06 alloys from first-principles 
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(Dated: July 25, 2011) 

We expand our study on cubic BiFeOa alloys presented in [K. Koumpouras and I. Galanakis, J. 
Magn. Magn. Mater 323, 2328 (2011)] to include also the BiMnOa and Bi2MnFe06 alloys. For 
the latter we considered three different cases of distribution of the Fe-Mn atoms in the lattice and 
six possible magnetic configurations. We show that Fe and Mn atoms in all cases under study 
retain a large spin magnetic moment, the magnitude of which exceeds the 3 ^b- Their electronic 
and magnetic properties are similar to the ones in the parent BiMnOs and BiFeOs compounds. 
Thus oxygen atoms which are the nearest-neighbors of Fe(Mn) atoms play a crucial role since they 
mediate the magnetic interactions between the transition metal atoms and screen any change in 
their environment. Finally, we study the effect of lattice contraction on the magnetic properties of 
Bi2MnFe06. 

PACS numbers: 75.50.Bb, 75.50.Ee, 75.70.Cn 



INTRODUCTION 

Spintronics brought to the center of the scientific re- 
search new materials with exotic properties [1|-|3|]. The 
latest addition to these materials are the so-called mul- 
tiferroics which combine several ferroic orders like fer- 
romagnetism, ferroelectricity, ferroelasticity etc [4|-l6|. 
Among them exist some compounds which combine elec- 
tric and magnetic order exhibiting the magnetoelectric ef- 
fect [4|. These alloys have several potential applications 
like magnetic-field sensors and electric-write magnetic- 
read random-access memories |7Hll|. Magnetic order and 
ferroelectricity have different origins |12h14 | and thus the 
materials exhibiting the magnetoelectric effect are few 
and the coupling between the magnetic and electric prop- 
erties is weak. An alternative route to achieve a strong 
coupling could be the growth of thin film heterostruc- 
tures and several advances have been made towards the 
magnetic control of ferroelectricity 
trie control of thin film magnetism 

Among the most studied single-component multifer- 
roic compounds are the bismuth ferrite and the bismuth 
manganite, BiFeOs and BiMnOa. Bulk BiFeOs crystal- 
lizes in a perovskitc-like pseudocubic structure instead of 
a ferrite one [2l| and bulk BiMnOs prefers a monoclinic 
lattice I22II. BiFeOs is a ferroelectric G-type antiferro- 
magnet |21| while BiMnOs is a ferromagnet presenting 
no net spontaneous polarization [23|, |24| . A Jahn- Teller 
distortion corresponding to an elongation of the oxygen 
octahedron surrounding the Mn atom induces a polar or- 
der in BiMnOa and G-type antiferromagnetism appears 
as in BiFeOs; this new multiferroic state is close in en- 
ergy to the ferromagnetic non-polar ground state |23l,l24|. 
Several first-principles calculations have been carried out 
to study the properties of both bulk BiFeOs [25| and 
BiMnOa compounds 23|, l2J, |26|. We refer readers to 
Ref. 27 for an overview of the literature on both com- 



present only a weak magnetoelectric effect, an alternative 
route to achieve a more strong effect has been proposed to 
be the growth of heterostructures where epitaxial strain 
can enhance the phenomenon [28-31.]. Doping has been 
also proposed to enhance the performance of such struc- 
tures 32ll33l|. Towards such heterostructures, multilayers 
consisting of alternating layers of BiFeOs and BiMnOa 
have been proposed where the ferroelectricity of BiFeOa 
couples to the ferromagnetic order in BiMnOa through 
epitaxial strain and a nanoscale checkboard from first- 
principles has been recently proposed for several lattice 
structures 2T\. Simple doping of BiFeOa with Mn leads 
to an antiferrornagnetic coupling of the Mn spin moments 



to the Fe ones [34 1 




and the elec- 



In a recent publication (Ref. |35I ) we have pre- 
sented extended first-principles calculations, employ- 
ing the Quantum-ESPRESSO i36| ab-initio electronic 
structure method in conjunction with the Generalized- 
Gradient Approximation (GGA) in the Perdew-Burke- 
Erzenhof formulation 371], on the electronic and mag- 
netic properties of BiFeOa alloy as a function of the lat- 
tice constant in the case of the cubic perovskite structure 
(see figure 1 in Ref. |35I ). where the heavy cations (Bi) 
occupy the corners of the cube, oxygen atoms the cen- 
ter of the faces and the Fe atom is at the center of the 
cube. Fe atoms alone form a cubic cell and the con- 
sidered 2x2 unit cell in our calculations contained eight 
primitive cells and thus eight Fe atoms. We studied four 
possible magnetic arrangements: one ferromagnetic (F- 
type) and three antiferromagnetic (A-type, C-type and 
G-type). We found that all four types of magnetic or- 
der are close in energy and for a lattice constant of the 
primitive cell larger than 3.888 A the G-type antiferro- 
magnetism becomes more stable than the ferromagnetic 
state. The density of states (DOS) and the spin mag- 
netic moments of the Fe atoms, which are responsible for 
the magnetic properties, showed similar behavior in all 
pounds. Since the single-component crystals like BiFeOa possible magnetic states since magnetic interactions are 
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TABLE I; Absolute values of the atom-resolved spin magnetic 
moments in jj,b for all cases and for all magnetic configura- 
tions under study. Results are for a lattice constant of the 
primitive unit cell (see Fig. 1 in Ref. [H) of 14.7 a.u. (3.888 
A). In case III the Al and A2 structures are equivalent due to 
symmetry reasons and the same is valid also for the CI and 
C2 configurations. 
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A2 

CI 

C2 
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BiFeOs BiMnOa 



I-case 



Bi2MnFe06 
Il-case 



Ill-case 



Fe 
2.98 
3.26 

3.53 

3.65 



Mn 
3.63 
3.53 

3.41 

3.24 



Fe 
2.99 
3.24 
3.26 
3.57 
3.51 
3.68 



Mn 
3.48 
3.49 
3.39 
3.32 
3.40 
3.22 



Fe 
3.18 
3.37 
3.30 
3.51 
3.58 
3.68 



Mn 
3.39 
3.27 
3.35 
3.37 
3.32 
3.21 



Fe 
3.31 
3.40 
3.40 
3.59 
3.59 
3.68 



Mn 
3.22 
3.26 
3.26 
3.31 
3.31 
3.18 



FIG. 1: Fe and Mn atoms in Bi2MnFe06 form a cubic lat- 
tice. In (a) we present the three distributions of Mn and 
Fe atoms which we considered and in (b) the possible mag- 
netic configurations (vectors pointing up correspond to pos- 
itive spin magnetic moments and vectors pointing down to 
negative moments). Notice that there are now two possible 
A-type antiferromagnetic structures: (i) Al where atoms be- 
longing to neighboring planes along the [001] direction have 
antiparallel spin magnetic moments, and (ii) A2 when this oc- 
curs along the [100] directions. For the same reason there are 
two C-type antiferromagnetic configurations along the [110] 
(CI) and along the [101] (C2) directions. 



mediated by oxygen atoms. Bi atoms showed no con- 
tribution to m.agnetism. while the magnetic properties of 
the oxygen atoms depended strongly on their local envi- 
ronment since they are located at the midpoints between 
neighboring iron atoms. 

We expand the study in Ref. ^ to cover also the case 
of Mn substitution for Fe in BiFeOa motivated by the 
work in Ref. |27| since the coexistence of Fe and Mn 

We 



cations seems to lead to new multifcrroic materials, 
have concentrated our interest on the electronic and mag- 
netic properties of the resulting compounds employing 
the same electronic structure method as in Ref. l35| and 
calculations details are similar. We have chosen to study 
a case with equal number of Fe and Mn atoms (denoted 
as Bi2MnFe06) using the same 2x2 unit cell as in Ref. 
35l . Thus now we have four Fe atoms and four Mn atoms 
per unit cell. There are several ways to distribute the 
Fe and Mn atoms and in panel (a) of Fig. [T]we present 
the three cases which we have considered. In the I-case 
we have a layered structure consisting of alternating pure 
Fe and pure Mn layers along the [001] direction. In the 
Il-case we have alternating pure Fe and pure Mn layers 
along the [110] direction and in Ill-case along the [111] 
direction. In I-case each Fe(Mn) atom has four Fe(Mn) 
and two Mn(Fe) atom as nearest transition-metal neigh- 
bors while in Il-case each Fe(Mn) atom has two Fe(Mn) 



and four Mn(Fe) atom as nearest transition-metal neigh- 
bors. In Ill-case each atom has six neighbors of the other 
chemical species. For our calculations we have used a lat- 
tice constant of 3.888 A for the primitive unit cell in Fig. 
1 of Ref. |35| which is the calculated equilibrium lattice 
constant for BiFeOa in [SSl Moreover due to the lower 
symmetry there is for Bi2MnFe06 a larger number of 
possible magnetic configurations with respect to BiFcOa 
which are presented in panel (b) of Fig. [TJ Arrows indi- 
cate the orientation of the spin magnetic moments of Fe 
and Mn atoms at the corners of the cube. F corresponds 
to the ferromagnetic configuration. There are two types 
of A-antiferromagnetism; in the Al-type all atoms within 
a (001) layer are coupled ferromagnetically but successive 
layers along the [001] direction are coupled antiferromag- 
netically between them while in the A2-type this concerns 
the layers along the [100] direction. CI- and C2-types of 
antiferromagnetism describe antiferromagnetic coupling 
of successive layers along the [110] and [101] directions, 
respectively. Finally G-type antiferromagnetism is char- 
acterized by alternation of the orientation of the spin 
magnetic moments along the [111] direction. We should 
note here that for the Ill-case of the distribution of Fe 
and Mn atoms, the symmetry is higher and the [001] and 
[100] directions are equivalent and the same stands also 
for the [110] and [101] orientations. Thus for this case 
Al and A2 as well as CI and C2 are degenerate. In the 
next section we present our calculated results. We have 
also calculated the case of BiMnOa similar to BiFeOa in 
Ref. |35| as reference and at the end of the next section 
we also present shortly results for a smaller lattice con- 
stant of 3.703 A corresponding to a contraction of the 
lattice of about 4.8 % . In the last section we conclude 
and summarize our results. 



RESULTS AND DISCUSSION 



BiFeO. 



BiMnO. 



We will start our discussion from the case of BiMnOa. 
In Table |I] we have gathered the spin magnetic moments 
in fiB for the transition metal atoms for all cases under 
study and for a lattice constant of the primitive cell of 
3.888A. We have also included in the first column the 
values for BiFeOa from Ref. [H Note that for BiMnOa 
and BiFeOs the Al and A2 as well as the CI and C2 an- 
tiferromagnetic configurations defined in Fig. [T] are de- 
generated. If we look at the spin magnetic moments Mn 
in BiMnOa posses in general a very high spin magnetic 
moment which reaches the 3.63 ^b in the ferromagnetic 
case and varies between 3.24 and 3.53 fiB in th^ antifer- 
romagnetic configurations comparable to the values for 
Fe in BiFeOa. The only noticeable difference for the two 
alloys concerns the ferromagnetic case (F-type) where Fe 
has a spin moment of about 3 fiB quite smaller than 
all other cases. Mn has one valence electron less than 
Fe and thus in order to understand the behavior of its 
spin magnetic moment with respect to Fe we have to 
look at the density of states (DOS) presented in Fig. [5] 
for Mn and Fe in BiMnOa and BiFeOa, respectively, for 
all four magnetic configurations under study. In the F- 
and A-configurations the large exchange splitting of the 
Mn d-states pushes the minority bands slightly higher 
in energy with respect to Fe while for the latter part of 
the extra electron occupies minority-spin states near the 
Fermi level leading to smaller Fe spin magnetic moments 
with respect to Mn. In the C- and G-configurations Mn 
presents a minority-spin pick pinned exactly at the Fermi 
level but this cannot explain the behavior of the spin 
magnetic moments in these cases. In the same figure we 
present with yellow shaded region the triple degenerated 
t2g states and thus the rest of the DOS corresponds to 
the double degenerated e^ states. The latter are higher 
in energy with respect to the t2g states due to the crys- 
tal field effect. The Cg orbitals (d22 and d2;2^j^2) point 
towards the neighboring oxygen atoms and thus ener- 
getically are unfavorable with respect to the t2g orbitals 
{dxy, dxz and dyz) which point in the intermediate space, 
since in the latter case the Coulomb repulsion with the 
electrons occupying the oxygen p-orbitals is smaller. In 
the minority-spin band of Mn the eg are not presented 
since they are slightly over the upper bound of the en- 
ergy axis. In the majority-spin band the eg states are 
partially occupied while for Fe the are completely occu- 
pied and thus Mn has a slightly smaller spin magnetic 
moment for the C- and G- configurations. Finally we 
should shortly discuss the other atoms in these two com- 
pounds. Bi atoms are heavy and their role is to yield the 
electric polarization when the lattice is distorted while 
due to symmetry their spin magnetic moment is almost 
zero and thus their contribution to the magnetic prop- 
erties can be neglected. Oxygen atoms as discussed in 
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FIG. 2: (Color online) For both BiFeOs and BlMnOa we 
present the Fe and Mn density of states (DOS) for all four 
possible magnetic configurations. With thick lines is the DOS 
projected on the d-orbitals and with the shaded region the 
contribution from the triple-degenerated t2g states. The re- 
maining contribution comes from the double-degenerated e^ 
states. Positive DOS values concern the majority-spin elec- 
trons and negative DOS values the minority-spin electrons. 
The zero energy value corresponds to the Fermi level. 



Ref. 35 for BiFeOa are located at the midpoints between 
neighboring transition-metal atoms. If the two Fe atoms 
had antiparallel moments the oxygen atom in the mid- 
dle has a zero spin magnetic moment due to symmetry, 
otherwise its spin moment was about 0.15 fiB due to the 
hybridization between the 0-p and Fe-t2g orbitals. Man- 
ganese t2g-orbitals hybridize with the oxygen p-orbitals 
much weaker and as a result oxygen atoms in BiMnOa, 
when both Mn neighboring atoms have parallel spin mo- 
ments, exhibit spin magnetic moments of one order of 
magnitude smaller than in BiFeOa (the larger obtained 
value is about 0.03 /is in the ferromagnetic configura- 
tion). Thus in BiMnOa the magnetic properties are en- 
tirely localized at the Mn sites. 

Probably the most interesting from the three cases of 
possible atomic distributions of Fe and Mn atoms pre- 
sented in Fig. [T] is the I-case where we have along the 



Bi^MnFeO^ (I-case) 
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FIG. 3: (Color online) Fe and Mn d-resolved DOS for the 
I-case (see atoms distribution if Fig. [T| for all six magnetic 
configurations. Notice that in all cases the positive DOS val- 
ues correspond to the majority states (e.g. in the Al and A2 
cases the majority states are the spin-up for the Fe atoms and 
the spin-down for the Mn atoms). 
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FIG. 4: (Color online) Same as Fig. [3] for Fe (upper panel) 
and Mn (lower panel) in Bi2MnFe06 for the II and III cases 
of atoms distribution if Fig. [l] 



[001] directions successive layers of pure Fe and pure Mn 
atoms. In Table U we have gathered the spin magnetic 
moments for all magnetic configurations under study and 
in Fig. |3] the Fe- and Mn- resolved DOS. Note in each 
magnetic configuration all Fe(Mn) atoms have the same 
magnitude of the spin magnetic moment (only its orien- 
tation varies) and thus the same DOS shape. We present 
the absolute value of the spin magnetic moments and 
in Fig. [3] the positive DOS values correspond to the 
majority-spin electrons in order to make comparison be- 
tween the various cases easier to follow. We will start our 
discussion from the spin magnetic moments presented 
in Table HI Our first remark is that both Fe and Mn 
spin moments are very close to their values in the par- 
ent BiFeOa and BiMnOa alloys, respectively. Moreover 
the Al and A2 as well as the CI and C2 antiferromag- 
netic configurations correspond to almost identical spin 
magnetic moments. Thus, as discussed in Rcf. 35,, we 
can conclude that the oxygen atoms play a crucial role. 
Although themselves posses very small spin magnetic mo- 
ments they bridge the transition-metal atoms and they 
mediate the magnetic interactions playing a shielding role 
for the Fe and Mn atoms. Each Fe(Mn) atom has now 4 



Fe(Mn) atoms as nearest transition-metal neighbors in- 
stead of 6 in the perfect alloys and the other two have 
been substituted by Mn(Fe) atoms. But this change in 
the environment is screened by the intermediate oxy- 
gen atoms and the spin moments of the transition-metal 
atoms are only marginally affected. This discussion ia 
also reflected in the DOS presented in Fig. [S] For each 
magnetic configuration the Fe- and Mn-resolved DOS are 
almost identical to the ones for the parent BiFeOs and 
BiMnOs alloys in Fig. [2] and all the details shown in Fig. 
[5] are present also in the DOS in Fig. |31 e.g. in both 
figures Fe in the F-case presents a double pick struc- 
ture in the minority-spin band and the Fermi level is 
pinned exactly between these two picks while in the G- 
configuration there two distinct minority-spin picks cor- 
responding to the t2g and eg states. 

For the II- and Ill-cases of the Fe-Mn distribution we 
include the spin magnetic moments in Table |T] and the 
DOS in Fig. IH We do not present the A2 and C2 cases 
in the figure since for the Il-case they are similar to the 
Al and CI magnetic configurations while for the Ill-case 
as mentioned in the introduction we cannot distinguish 
between the Al(Cl) and A2(C2) cases due to the sym- 
metry of the lattice presented in Fig. [TJ In Il-case we 



have [001] chains of pure Fe and pure Mn or equivalently 
we can describe the structure as successive layers of pure 
Fe and pure Mn along the [110] direction. With respect 
to the I-case now each Fe(Mn) atom has only 2 nearest- 
neighboring transition metal atoms of the same chemical 
species instead of 4. In the Ill-case, which can be en- 
visaged as alternating layers of pure Fe and pure Mn 
along the [111] direction all nearest neighbors are of the 
other chemical type, e.g. each Fe atom has 6 nearest- 
neighboring Mn atoms. Thus although the magnetic in- 
teractions are mediated by the oxygen atoms the effect 
on the spin magnetic moments is larger especially for the 
ferromagnetic (F-type) alignment of the spin moments. 
As the number of Mn neighbors increases the Fe spin 
magnetic moment in the F-case increases from 2.99 /is 
in the I-case to 3.18 /is in the Il-case and 3.31 /is in 
the Illl-case. Simultaneously the Mn spin magnetic mo- 
ment drops from 3.48 jjlb in the I-case to 3.22 ^b in the 
Ill-case. For the antiferromagnetic cases under study the 
influence on the spin magnetic moments is smaller and for 
the G-type antiferromagnetism spin magnetic moments 
are unaltered by the variation of the Fe-Mn distribution. 
Overall even in the F-case the change in the magnitude 
of the spin magnetic moments is quite small; less than 0.3 
^B ■ The behavior of the spin magnetic moments is also 
reflected on the DOS presented in Fig. 21 In the Il-case 
both Fe- and Mn-resolved DOS for all magnetic configu- 
rations under study are similar to the I-case. In the III 
case the only noticeable difference concerns Fe in the fer- 
romagnetic case where the two picks at the Fermi level in 
the minority-spin band have merged in one broad large 
pick due to the influence of the exclusively Mn neighbors 
and this is also the only case where the variation of the 
spin-magnetic moment is noticeable reaching the 0.3 /is 
with respect to the I-case. 

Finally, we have also studied the effect of a smaller 
lattice constant for the I-case and in Table HI] we present 
the obtained results for a lattice constant of 3.703 A, 
which is about 4.8 % smaller than the 3.888 A used in all 
previously presented results. Already in Ref. |35| we had 
shown that contracting the lattice of BiFeOs leads to an 
important decrease of the spin magnetic moment which 
in the G-type antiferromagnetism is less than 1 /ib with 
respect to the 3.65 hb in Table U for the larger lattice 
constant. Also in the case of Mn in BiMnOa the Mn 
spin magnetic moment decreases with the contraction of 
the lattice although the decrease is not as dramatic as 
in BiFeOa. This picture is also reflected in the layered 
I-case of Bi2MnFe06 which we studied. In the I-case the 
Fe atoms keep a large number of Fe nearest-neighbors 
and its spin magnetic moment exhibits an even larger 
decrease than in the case of BiFeOa. Especially in the G- 
type antifferomagnetism the Fe spin magnetic moments 
attend a low of only 0.25 ^b being at the edge of being 
non-magnetic. On the contrary Mn atoms keep a high 
value of their spin magnetic moment which exceeds the 



TABLE IL Same as Table|T]for a lattice constant of the primi- 
tive cell of 14 a.u. (3.703 A) corresponding to a contraction of 
the lattice constant of about 4.8 % with respect to the results 
in Table m 
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2.32 


3.05 


2.21 


3.08 
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2.01 


3.03 


1.70 
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1.87 


3.05 


CI 


1.55 


2.93 


1.07 


2.81 


C2 






0.76 


3.11 


G 


0.98 


2.53 


0.25 


2.84 



2.8 /is in all cases and thus, as in the case of BiMnOa, Mn 
magnetic properties are less affected by the contraction 
of the lattice. 



SUMMARY AND CONCLUSIONS 

We expand our study in Ref. i35 to the case of 
Bi2MnFe06 alloys. First we have studied the BiMnOa 
alloy for various magnetic configurations and we have 
shown that, as was the case for BiFeOa, Mn atoms exhibit 
high values of their spin magnetic moment, the magni- 
tude of which exceeds the 3 /is in all cases under study. 
The small difference between the Fe and Mn spin mag- 
netic moments can be easily explained by the position of 
the eg states with respect to the t2g states which lie lower 
in energy due to the crystal field effect. We took into ac- 
count three different possible distributions for the Fe and 
Mn atoms in Bi2MnFe06 and studied one ferromagnetic 
and five possible antiferromagnetic alignments of the spin 
magnetic moments. Both Fc- and Mn-resolved magnetic 
and electronic properties are only marginal sensitive to 
the various distributions of the atoms and show prop- 
erties close to Fe and Mn atoms in the parent BiFeOa 
and BiMnOa alloys, respectively. Also among the vari- 
ous magnetic configurations the spin moments show small 
variation and their magnitude is always larger than 3 
fiB ■ Thus oxygen atoms which are the nearest- neighbors 
of Fe(Mn) atoms play a crucial role since they mediate 
the magnetic interactions between the transition metal 
atoms and they screen any change in their environment. 
Finally we have shown that contracting the lattice leads 
to a large decrease in the spin magnetic moments of the 
Fe atoms as in BiFeOa studied in Ref. l35| while the Mn 
atoms retain a large portion of their spin magnetic mo- 
ment. 
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